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Abstract 

By using a string-inspired modular invariant supergravity model, which was proved well 
to explain WMAP observations appropriately, a mechanism of preheating just after the end 
of inflation is investigated. By using the canonically normalized and diagonalized scalars, 
the decay rates of these fields are calculated inflaton S into gauge sector fields. 

The reheating temperature is estimated by both the stability condition of Boltzmann 
equation and the instant preheating mechanism. Both of reheating temperatures are almost 
the same order of magnitude ~ O(10 10 ) GeV. Because two mechanisms are completely 
independent processes, the former is caused through the inflaton decays into gauge fields 
and gauginos and the latter is caused by the scattering process of two inflatons into two right 
handed sneutrinos, which will decay into Higgs fields and other minimal SUSY standard 
model (MSSM) particles, we conclude that both mechanisms play essential roles in the 
preheating process after inflation. 



1 Introduction 



Following "Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations" [JJ, the 
theory of inflation are proved to be the most promising theory of the early universe before the 
big bang. 

As a favored scenario to explain the observational data, it is customary to introduce a scalar 
field called inflaton. What kind of theoretical frameworks are the most appropriate to describe 
the theory of particle physics, inflation and the recent observations? It seems to require a 
far richer structure of contents than that of the standard theory of particles. As far as the 
AD, N = 1 supergravity can play an elementary role in the theory of the space-time and the 
particles [5] , it can also be essential in the theory of the early universe as an effective field theory. 
Supergravity, however, has been confronting with the difficulties, such as the 77-problem and the 
supersymmetry breaking (SSB) mechanism has been studied by many authors [2 |H 02 [6] . We 
have investigated to prevail over these difficulties in Refs.[71|5] by using the modular invariant 
supergravity induced from superstring [5]. We found that the interplay between the dilaton 
field S and gauge-singlet scalar Y could give rise to sufficient inflation. The model is free from 
the the 77-problem and realizes appropriate amount of inflation as well as the TT angular power 
spectrum. 

In this letter, the preheating mechanism just after the end of inflation will be investigated. 

First we will briefly review the model and the former results [TJIEj as follows. It is convenient 
to introduce the dilaton field S, a chiral superfield Y and the modular field T. Here, all the 
matter fields are set to zero for simplicity. Then, the effective No-Scale type Kahler potential 
and the effective superpotential that incorporate modular invariance are given by [9], (see also 

my 

K = -hi(S + S*)-3 In (T + T* - \Y\ 2 ) , (1) 
W = 3bY 3 In [c e s/3b Yrj 2 (T)l , (2) 

where rj(T) is the Dedekind 77-function, defined by: 

00 

7?(T) - e -^/2i -Q (1 _ e -2,n T) (3) 
n=l 

The parameter b and c are treated as free parameters in this letter as discussed in Ref.[5]. We 
had found that the potential V(S, Y) at T = 1 has a stable minimum at for the values b — 9.4, 
c = 131 and obtained 

S min = 1.51, r min = 0.00878480, (4) 

where r)(l) = 0.768225, ?y 2 (l) = 0.590170, r)'(l) = -0.192056, 17" (1) = -0.00925929 are used. 
The inflationary trajectory can be well approximated by 

Y min (S) ~ O.009268e-° 035461S , (5) 

which corresponds to the trajectory of the stable minimum for both S and Y. The slow- 
roll parameters es and rjss satisfy the slow-roll conditions. The number of e-folds ~ 57, by 
integrating from S cn d ~ 4.19 to ~ 11.6, i.e. our potential can produce a cosmologically 
plausible number of e-folds lj. Here is the value corresponding to the scale &;* = 0.05 
Mpc -1 . We can also compute scalar spectral index and its running that describe the scale 
dependence of the spectrum of primordial density perturbation Vn = (H/S) 2 (H/2tt) 2 [TT] : 
these indices are defined by 

ns - l = ^k> (6) 
a = dn s 
din A: 
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These are approximated in the slow-roll paradigm as 



n a (S)~l-6e s + 27755, (8) 
a s (S) ~ 16e sVss - 24e| - 2£ ( 2 3) , (9) 

where is an extra slow-roll parameter that includes the trivial third derivative of the poten- 
tial. Substituting 5* into these equations, we have rt s * ~ 0.951 and a s * ~ —2.50 x 10~ 4 . 

Because n s is not equal to 1 and a s is almost negligible, our model suggests a tilted power law 
spectrum. The value of n s » is consistent with the recent observations; the best fit of seven-year 
WMAP data using the power law ACDM model is n s ~ 0.963 ± 0.014 [1 . Finally, estimating 
the spectrum V-r. in slow-roll approximation (SRA), 

^l^rW' (10) 

we find V-r* ~ 2.36 x 10 -9 . This result matches the measurements as well [T][7J[5]. Incidentally 
speaking, the energy scale V ~ 10~ 10 GeV is also within the constrained range obtained by 
Liddle and Lyth [TJ. 

In order to study on the angular power spectrum, we need the tensor perturbation (the 
gravitational wave) spectrum which is given as follows: 

H \ 2 2 



In SRA, the spectral index of "P gra v is given by the slow-roll parameters e and 77 as 

n T = -2e. (12) 
Using these parameters TT and TE angular spectrum were well fitted to the WMAP data [3]. 

2 Gravitino mass and the other mass parameters 

Now we will briefly investigate the properties of inflaton S, gravitino and SSB mechanism. First, 
gravitino mass is given in this case 

m 3/2 = M P e K/2 \W\ = 3.16 x 10 12 GeV, (13) 

where h = 6.58211915 x 10~ 25 GeV-sec and M p = 2.435327 x 10 18 GeV are used. 

In our model, by expanding the potential V around the minimum of S(t), Y(t) and fixed 
T = 1, and by providing S(t) and Y(t) are real, then we obtained S(t), Y(t) as follows: 

S(t) = S mia + ^-^l, (14) 
V o mst 

! '") = ?(i^ e "*' (15) 

After scalars S, Y, T are canonically normalized and the masses diagonalized [13] , [14] , the mass 
eigenstates are denoted by S',Y',T', then masses are calculated as Mgi = 3.97 x 10 12 GeV, 
My/ = 2.45 x 10 17 GeV, M T > = 9.02 x 10 12 GeV, where S', Y', T are defined as follows: 

5" = 3.00 x lO^S" + 1.94 x 10~ 3 y - 3.66 x lO" 1 ! 1 (16) 

Y' = 3.82 x lO" 4 ^ + 1.22 Y - 2.94 x 10" 8 T (17) 

T' = 1.40 x lO" 1 ^ - 7.49 x 10" 3 V + 7.85 x lO^T. (18) 
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Supersymmetry is overwhelmingly broken by superficial S", which will be shown in separate 
paper comparing with the fact pointed out by Nilles et al. jT5J [TH1 [T7] , in which the interchange 
of SSB fields occurs (see also Kalolosh et. al. [TB]). 

Canonically normalized fcrmionic states of supersymmetric partners S,Y,T are given by 

S" = 0.331S", Y' = 1.22V, f ' = 0.867T - 7.61 x 10~ 3 Y, (19) 

and the values of them are numerically determined as 

mg, = 0GeV, m y , = 3.01 x 10 17 GeV, m f , = 2.65 x 10 15 GeV. (20) 

Since S is massless and S breaks supersymmetry, S state is identified with Goldstino, which is 
absorbed into gravitino by super-Higgs mechanism [3J [TS] . 

Non-thermal production of gravitinos is not generated from the inflaton (dilaton), since the 
inflaton mass is lighter than gravitino, but they are produced by the decay of modular field T 
and scalar field Y , which processes are shown in our separate paper. 



3 Decay rate from inflaton to gauginos 

Because the canonically normalized mass eigen state inflaton 5" does not decay into gravitinos, 
S' will decay directly decay into the minimal SUSY standard model (MSSM) particle or the 
next to minimal SUSY standard model (NMSSM) particles after the end of inflation. As an 
example, the decay rate of S' into gauginos is estimated in our model. By using the term 
Cgaugino = «/ d 2 0f a b(6)W a W a , f a b{4>) = <fr5 a b, the interaction between S and gauginos A a 's 



are given as 
r 



2/3(0) 



l dfab{4>) K/2, 

~4 86 



G lPr D r W*X a X b 



-fi b (6)V, [AVA 1 ] 



e K l 2 G^ r D^WX a X b . (21) 



By seeing the first term of (|2"Tj) , A°'s are also canonically normalized as A" = (f^ b ) 2 A ■ 

The interactions come from the third and fourth terms. The terms include e K / 2 G^ D^W* 

which implies the auxiliary field of <j) in global SUSY theory and it is replace by F^. 

By expanding ^g^-F^ in the terms around the stable point, interaction terms are given as 



£int — 



4(/ab) 



d 2 fat 



F„ 



dfab dFj, 
86 86 



66- 



dfab 8F^ 



4 (fab) L 



d 2 rab F ^df; lb 8F;, 



86 86 

dflb dF r 



a \ b 



x a x 



(22) 



where when 6 = S, F$ implies the SSB scale of the model and will be estimated as (S + S*} 
to 3 /2 since (Fg) ~ m 3 / 2 and (S+S*) take value about 3 times of Planck scale. Therefore, as the 



first term contribute far smaller than the second and negligible, — (^&) 



1 3/2 



is remained. 



The derivative term by S* can be replaced by 
A + A) can be estimated as: 



(aP") ~ ms - Then the decay rate T( 



T(S -)• A + A) 



(«i> 2 



t 3/2 



167r (fab) 



, "^3/2 

' m s 



4m\ 

771% 



(23) 



By using the relation Fg ~ M p mgp that holds for the mass of SUSY particles, the order of 
gaugmo mass m\ will be estimated as O(10 7 ) GeV - O(10 8 ) GeV. Then the decay rate of 
inflaton to gauginos can be estimated to be T(S — > X + A) ~ 3.89 x 10 3 GeV [2U] . 
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4 Reheating temperature 



Now we will calculate the reheating temperature from two methods and compare the resulting 
temperatures. The first method is using Boltzmann equation and the decay rate of the inflaton 
already calculated above. Another rely on the instant parametric resonance to calculate the 
evolution of the number density of the produced matters and obtain reheating temperature. 

First, the reheating temperature Tfj(gaugino) is derived from Boltzmann equation by using 
the decay rate, is given by 

10 



^(gaugino) = y/Mp T(S — > A + A), (24) 

and numerically given as T R ~ 4.45 x 10 10 GeV, by inserting the decay rate from the canonically 
normalized inflaton field S'. where g* is the number of the effective degrees of freedom of MSSM, 
i.e. = 228.75. 

Next, we will estimate the reheating temperature T R by instant parametric resonance by in- 
troducing NMSSM superpotential with the righthanded neutrino superfields. We should choose 
a model to determine the reheating temperature by inflaton decay into MSSM or NMSSM 
particles. We will assume the instant preheating mechanism |21j . because this method is math- 
ematically easier to control than that of parametric resonance |22) . 

It is not unique to take in the ordinary particles into string-inspired modular invariant 
supergravity. We will assume that minimal Kahler potential J^i ^t®* 1 is simply added to K 
and super potential of NMSSM is added to include a term that directly coupled with the inflaton 
(dilaton) superfield during the oscillations of inflaton. We will choose it as that the righthanded 
neutrino superfields N c couple directly with the inflaton [25] : 

Wnmssm = M Ri N?N? + \iSN9N? + rfj N?LjH u . (25) 

Number density through the parametric resonance is given by |22j 



4 +1 



"™ 2 + (l + 2e"™ 2 ) n{ - 2e~™ 2/ Vl + e-™ 2 \Jn{(l + n{) sin 9 j . (26) 



If we consider the instant preheating mechanism |24j , the number density is simply given by the 
first term 

7,(k 2 /a 2 + M R . ) 

n k = e-™ 2 = e ^ , (27) 

where we have used NMSSM superpotential (|25[) . By integrating this equation in k, we derive 
the number density of supersymmetric partner of the righthanded neutrinos n ^ a . 

T2 



~ 2,* J Q dkknk = ^^ eXP {-^' (28) 



From this equation we can estimate the reheating temperature 

f2 



1/4 



T R (l^nt) = [ -s- • M Ri • n* f j ^ exp j . (29) 



Since = 6.39 x 10 in our model, only free parameters in T^(Instant) are Ai and M Ri . We 
assume here to restrict Xi — Af^./|So| so as to the index of the exponent in the number density 



Ri 

takes the value O(l), then 



1/4 



T fl (Instant) = M Ri ( — =— e" 77 I , (30) 
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Tr (Instant) only depends on M^. If we estimate the value of Mr ( ~ ms as O(10 12 ) GeV, 
takes the value O(10 -6 ) and Tfj(Instant) becomes O(10 10 ) GeV, which value is similar value as 
Tfl(gaugino). Therefore we conclude that both the contribution from inflaton decays and that 
from the parametric resonance of four body scattering process play equally important roles in 
the preheating process just after the end inflation. 

Because the primordial gravitinos decay very rapidly and the reheating temperature is lower 
than the gravitino mass, the effect to the standard Big Bang Nucleosynthesis (BBN) scenario 
[U [53 W\ M W\ ISO]) may be negligible in our model (see also [23 131])- 

5 Conclusion 

We have investigated on the preheating mechanism just after the end of inflation through both 
the inflaton (dilaton) decay into MSSM gauge sector and the collision of two inflaton into two 
righthanded sneutrinos. We conclude that the contribution of both the inflaton decays and 
the parametric resonance of four body scattering process play equally important roles in the 
preheating process just after the end inflation. The model we used, cleared the 77-problem and 
appeared to predict successfully the values of observations at inflation era. It predicted for 
examples, the indices n s>t ~ 0.951 and a s » ~ —2.50 x 10~ 4 . The value of n sie is consistent with 
the recent observations; the best fit of seven-year WMAP data using the power law ACDM 
model is n s ~ 0.963 ± 0.014 [1 . The estimation of the spectrum was as Vk* ~ 2.36 x 10~ 9 , 
which result matches the measurements as well [TJ [8] . 

The reheating temperature Tr is estimated by assuming the instant preheating mechanism 
and by using a tentative model among NMSSM models and comparatively the same order 
with that from the decay process of inflaton into MSSM gauge sector, which value is about order 
~ O(10 10 ) GeV. Because the mass of gravitino is calculated as 3.16 x 10 12 GeV, it is rather heavy 
and may be unstable, therefore, may not be considered as the lightest supersymmetric particle 
(LSP) or the next lightest (NLSP) and not a dark matter candidate discussed in Refs. [13lH4ll33] . 
However the main topic of supergravity at present stage of the theory is whether the gravitino 
exist or not in nature despite its mass. It is not reproduced after the reheating of the universe. 
The gravitinos are produced almost instantly just after the end of inflation through Y and T, 
not from inflaton. Then the yield variable for gravitino may take rather large value, however 
the decay time appears very rapid and disappear before the BBN stage of the universe. Though 
the effects of this type of gravitinos in the evolution of the universe should be investigated more 
carefully, the topic must be remained to later works. Therefore, we only remark here that our 
present model seems consistent with the present situation of observations. 

On the other hand, we commented that supersymmetry is overwhelmingly broken by F— term 
of the inflaton (dilaton) superfield S, that may be contrary to the occurrence of the interchange 
of SSB fields pointed out by Nilles et al. [T31 [T5J |T7] . More detailed investments will be shown 
in our later works. 

Though we have been exclusively restricted our attention to a model of Ref.[S], the other 
models derived from the other type of compactification seems very interesting. Among them 
KKLT model [3H [351 OS [37] attracts our interest, where the moduli superfield T plays an 
essential roles. We should take all the circumstances into consideration on essential problems 
confronted in construction of (string-inspired modular invariant) Supergravity models. 
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